Microcanonical phase transitions in small systems 
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When studying the thermodynamic properties of mesoscopic systems the most appropriate mi- 
crocanonical entropy is the volume entropy, i.e. the logarithm of the volume of phase space enclosed 
by the hypersurface of constant energy. For systems with broken ergodicity, the volume entropy has 
discontinuous jumps at values of energy that correspond to separatrix trajectories. Simultaneously 
there is a convex intruder in the entropy function and a region of negative specific heat below such 
critical energies. We illustrate this with a simple model composed of a chain of 3 particles which 
interact via a Lennard- Jones potential. 
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I. INTRODUCTION 

The work of D.H.E. Gross [l[ has recently pointed out 
the attention on the fact that a microcanonical descrip- 
tion of systems which may display phase transitions is 
in general more adequate than the traditional canoni- 
cal one. This is because the canonical description may 
"smear out" important information contained in the mi- 
crocanonical description which is richer [l[. For exam- 
ple negative specific heats, which have recently been ob- 
served experimentally in mesoscopic systems 0, H, H[, 
can be accounted for in the microcanonical ensemble but 
not in the canonical one Indeed it is well known 
that canonical ensemble and microcanonical ensemble are 
not in general equivalent, even when the thermodynamic 
limit is considered [J, @, 0, [|[ . 

The statistical mechanical analysis of physical sys- 
tems based on canonical ensemble is quite well estab- 
lished and universally agreed upon. Roughly speaking 
one has to compute the partition function Z{(3) and de- 
rive the thermodynamics of the system from the free en- 
ergy F = — /3 _1 In Z([3). Things are not quite so broadly 
agreed upon in the case of the microcanonical ensem- 
ble. In fact since the pioneering works of Boltzmann and 
Gibbs two possibilities were given for the microcanon- 
ical analysis of physical systems, which correspond to 
two different definitions of entropy (see for example the 
textbook of Gibbs [H or the more recent textbook of 
Huang Following [ll[ and [12] we shall refer to 

these two entropies as "surface entropy" and "volume 
entropy". The surface entropy is defined as: [3^ 



where 



Sa(E) = lnQ(E) 



Q(E) = J dzS(E-H(z)) 



(1) 



(2) 



with 5(x) denoting Dirac delta function. Sometimes this 



is also referred to as Boltzmann entropy or Boltzmann- 
Planck entropy. The volume entropy is defined as: 



where 



S*(E) = ln$(£) 



= J dz9(E - H{z)) 



(3) 



(4) 



'Electronic address: campisi@unt.edu 



with 6(x) = dyS(y) denoting Heaviside step func- 
tion. The quantities f2 and <& are related through the 
differential equation [l3j |: 

$'(£) = n(E) (5) 

where the prime symbol denotes derivation with respect 
to E. These entropies are named surface entropy and 
volume entropy because they are calculated as the loga- 
rithm of the area of the hyper-surface of constant energy 
in phase space and the volume of phase space that it 
encloses respectively. 

With reference to the literature about microcanonical 
phase transitions the surface entropy is certainly the most 
popular. For example Barre et. al. 0] method based 
on large deviation techniques uses the surface entropy. 
The surface entropy is used also in Rugh's microcanoni- 
cal formalism adopted in Ref. The strongest advo- 
cate of surface entropy is perhaps Gross [HI]. Nonethe- 
less pioneers of microcanonical phase transitions, such as 
Thirring [lfl ] and Lynden-Bell p|, used the volume en- 
tropy. 

The two entropies coincide in the thermodynamic limit 
but when the number of degrees of freedom of the sys- 
tem under study is small relevant differences may appear, 
therefore it is necessary to choose properly. Some au- 
thors [13, G3] have already pointed out that the surface 
entropy is not adequate when dealing with small systems 
because it does not account properly for finite-size ef- 
fects. On the other hand there is a number of theoretical 
reasons to prefer the volume entropy when the number of 
degrees of freedom becomes small. Here we shall review 
these reasons and we will illustrate the employment of 
the volume entropy in the study of phase transition with 
a small Lennard- Jones chain which displays a region of 
negative heat capacity. 
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II. WHY VOLUME ENTROPY 

In this section we will summarize some old and recent 
results concerning the volume entropy. These results in- 
dicate that, no matter the number of degrees of freedom 
of the system under study, the volume entropy always 
provides a good mechanical analogue of thermodynamic 
entropy. For example Helmholtz [l8[ proved that the 
logarithm of the area enclosed in phase space by the tra- 
jectory of a 1-dimensional system (i.e., the ID volume 
entropy «S$) provides a mechanical analogue of physical 
entropy in the sense that if one considers the quantities 



P 



/ dH \ 



and T = (2K) t where K is the kinetic energy, 
V is an external parameter on which the Hamiltonian de- 
pends and (-) t denotes time average, then, 



dE + PdV 
T 



exact differential = dS$ 



(6) 



This result, known as Helmholtz Theorem 0, [uj, says 
that the volume entropy is a good mechanical analogue 
of thermodynamic entropy in the sense that it repro- 
duces exactly the fundamental law of thermodynamics 
(i.e. the heat theorem ©)■ The Helmholtz Theorem 
has been recently generalized to multi-dimensional er- 
godic (i.e. metrically indecomposable) systems (see Ref. 
flit or Ref. [l!| for a different but equivalent approach). 
The resulting Generalized Helmholtz Theorem essentially 
states that the volume entropy reproduces the heat theo- 
rem no matter the number of degrees of freedom N. The 
same cannot be said about surface entropy which has 
been proved to reproduce it only up to corrections of the 
order 0(1/N) @. As a matter of fact Gibbs presents the 
volume entropy in his celebrated Principles of Statisti- 
cal Mechanics Q as the entropy that naturally satisfies 
the fundamental principle of thermodynamics (that is the 
heat theorem) |20j. 

Hertz pl| pointed out that the volume entropy is 
an adiabatic invariant already in 1910, and based his 
approach to statistical mechanics on it. Among the 
textbooks that adopt the same approach, those of 
Minister [12], Becker (23|, and the more recent book of 
Berdichevsky [24j are worth mentioning. Adiabatic in- 
variance is another good property of volume entropy be- 
cause it reproduces quite well Clausius' requirement that 
"For every quasi static process in a thermally isolated 
system which begins and ends in an equilibrium state, 
the entropy of the final state is equal to that of the ini- 
tial state" [25]. Of course the surface entropy is not an 
adiabatic invariant, although it becomes approximately 



such as the number of degrees of freedom increases [11 



Very recently it has been also proved that non- 
adiabatic transformations occurring in isolated systems 
which arc initially in a state of thermal equilibrium al- 
ways result in an increase of the expectation value of the 
volume entropy [26|, [27[. This result too does hold no 
matter the number of degrees of freedom N and cannot 
be proved in general for surface entropy. Thus the volume 
entropy explains quite satisfactorily also Clausius' law of 



entropy increase "For every non quasi static process in 
a thermally isolated system which begins and ends in an 
equilibrium state, the entropy of the final state is greater 
than that of the initial state" [25| . 

Recently more and more authors are becoming aware 
of the theoretical value of volume entropy. For example, 
on the basis of a Laplace transform technique for the mi- 
crocanonical ensemble, Pearson et. al. [281 ] reached the 
conclusion that the volume entropy "is the most correct 
definition for the entropy, even though it is unimportant 
for any explicit numerical calculation" , meaning that in 
the thermodynamic limit the difference with surface en- 
tropy becomes negligible. On the other hand for small 
system, such intrinsic correctness of the volume entropy 
becomes very important. Adib [l2j argues that the fi- 
nite size corrections to surface entropy found in Ref. (29| 
would be unnceded if the volume entropy were used in- 
stead. 

It is worth mentioning that the volume entropy has 
another property that is particularly important for small 
systems which have negative heat capacity, namely it is 
a naturally nonextensive entropy. According to Lynden- 
Bell Q, in fact, systems with negative heat capacity are 
necessarily nonextensive. The property of nonextensivity 
of volume entropy follows directly from the composition 
rule of enclosed volumes, $i i = 1, 2 , of two systems with 
total energy E = E x + E 2 [3: 



$(E)= / dEa$i(£a)$ 2 (£-£aJ 



(7) 



which is not a simple multiplication but a form of convo- 
lution which accounts for all possible partitions of ener- 
gies between the two systems. It has to be stressed that, 
despite of what is often stated in literature, the surface 
entropy is nonextensive too, as the composition rule for 
surface integrals Qi is the convolution, not the multipli- 
cation IFJ: 



n(E)= I 'dE 1 ni(JSi)n 2 (JS - Et 



(8) 



In sum, the volume entropy accounts for certain basic 
principles of thermodynamics, like the heat theorem and 
Clausius formulation of the second law equally well for 
large and small systems, whereas the surface entropy ac- 
counts for them only in the case of large system. For this 
reason it is the most appropriate mechanical analogue 
of thermodynamic entropy when dealing with small sys- 
tems. 



III. LENNARD-JONES CHAIN 

According to the Helmholtz Theorem [llj the mechan- 
ical analogue of thermodynamic entropy of a one dimen- 
sional system is 



S*{E, V) = log2 / dxy / 2m (E-(p(x,V)) (9) 

Jx-(E,V) 
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FIG. 1: Phase space structure for a particle of mass m — 1 
in a Lennard- Jones box of size V = 4 > V c . The separatrix 
corresponds to the critical energy E c = —0.0308. Below the 
critical energy there are two distinct trajectories (the dynam- 
ics is not ergodic). 
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FIG. 2: Entropy versus Energy for a particle of mass m = 1 
in a Lennard- Jones box of size V = 4 > V c . The discontinuity 
of Entropy at the critical energy signals a discontinuous phase 
transition. 



where x± {E, V) denote the turning points of the tra- 
jectory. If the potential is such that there is only one 
trajectory per energy level (crgodicity), then S$ satis- 
fies Eq. ([6]) @, [Tl|. Nonetheless, if the system has more 
than one trajectory per energy level for a certain range 
of energies and the system is on one of them, still we 
can use the above formula and say that the heat theorem 
is satisfied as long as the energies considered lye within 
that energy range. In this case P and T would be calcu- 
lated as time averages over the actual trajectory and 5$ 
would be given by the area enclosed by that trajectory 
only. Let us illustrate this with a practical example. Let 
us consider a ID chain composed of three particles which 
interact via a Lennard- Jones potential. Let us fix the 
position of two of them and let us place the third one in 
between, so that the first two particles act as walls of a 
ID box. Let us now study the behavior of the particle 
inside the box. Let the interaction potential be: 



u{x) 



1 

Zl2 



(10) 



and let us place the "walls" at x = ±V/2. Then the 
particle in the box is subject to the following potential: 



<p(x, V) = u{x + V/2) + u(x - V/2) 



(11) 



For values of V larger than a certain critical value V c ~ 
2.5, this system has a critical energy E C (V) = 0(0, V) 
such that for energy below E C (V) ergodicity is broken 
and there are two trajectories per energy level. Above E c 
the dynamics is ergodic and there is only one trajectory 
per energy level. Figure [1] shows a contour plot of various 
energy levels in phase space for a particle of mass m — 1 
in the Lennard-Jones box of size V = 4 > V c . For 
energy E = E c we have a separatrix. Below E c the curve 
of constant energy splits into two disconnected curves, 
whereas for values of E larger than E c we have only one 
curve. Below the critical energy the volume integral $ 
is given by the area enclosed by one of the two possible 
trajectories. As the energy crosses the critical value the 




energy 

FIG. 3: Temperature versus Energy for a particle in of mass 
m = 1 in a Lennard-Jones box of size V = 4 > V c - The curve 
displays a region of negative specific heat. At the critical 
energy the temperature goes to zero. 



integral <!> jumps discontinuously. In formulae we have: 



-9(E C -E) + 6{E - E c 



x / [2m{E - ip{x,V))] 1 ' 2 dx (12) 



HE) 



1/2 

The symbol [y] , denotes a function that is equal to 
yjy for y > and is null otherwise. The volume entropy, 
which is calculated by taking the logarithm of the expres- 
sion above, then displays a jump at the critical energy as 
well. Figure [2] shows a plot of S as a function of E, for 
the values m = 1 and V = 4 > V c . The critical energy is 
E c = —0.0308. Figure [3] shows the temperature plotted 
against the energy calculated, according to Eq. ([6]), as 



T=(2K) t 



HE) 

n(E) 



dS* 
~dE 



(13) 



There is a region of negative slope in the graph which 
correspond to a negative heat capacity. 
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FIG. 4: Entropy versus Energy for E < E c . Right below the 
critical energy the entropy function (thick line) is convex (the 
thin straight line is only a guide for the eye). 



IV. DISCUSSION 

The example provided in the previous section is per- 
haps too simple to be of interest to any specific physi- 
cal problem. Nonetheless it illustrates qualitatively the 
mechanism of microcanonical phase transition as cap- 
tured by the volume entropy. Such phase transitions 
are associated with the crossing of separatrix trajecto- 
ries, for which the dynamics of the system has no finite 
time scale. The figures show neatly that at the separatrix 
energy the entropy has a discontinuous jump, the tem- 
perature goes to zero, and for energies below the critical 
energy we have a region of negative heat capacity. These 
are not specific features of the system studied [36|. When- 
ever a separatrix is crossed there is a sudden open-up of 
a larger portion of phase space for the trajectory to en- 
close which leads to a discontinuity in the entropy 37 1. 
Further, at the separatrix, the period of motion, which 
for a well known theorem of classical mechanics is given 



by $' = Q [3(j, becomes infinite. Therefore the temper- 
ature, i.e., T = goes to zero. Since the temperature 
T = 2 < K >t=< p 2 /m >t is definite positive, below 
the critical energy there necessarily is a region of nega- 
tive slope, that is negative heat capacity. The appearance 
of a negative heat capacity is associated with a convex 
intruder in the entropy (see Fig. [5]) which signals the 
approach to the separatrix from below. It is important 
to notice that using the surface entropy would lead to 
a drastically different result. In this case the tempera- 
ture would be calculated as Th = , which might not 
tend to zero at the critical energy! Note also that Th is 
not proportional to the average kinetic energy and can 
be negative. Therefore, in agreement with Ref. (l7| we 
believe that surface entropy is not suited for low dimen- 
sional systems with broken ergodicity. 

The volume entropy could be used to address micro- 
canonical phase transitions in small dimensional systems 
with either long or short range interactions, like the c^ 4 
model, chains of particles interacting via Lennard- Jones 
potential [3l[ or the Hamiltonian Mean Field model [12] • 
All these models are expected to undergo a breaking of 
ergodicity [jjH 133^ thus there are separatrix trajectories 
and possible phase transitions that the volume entropy 
can detect. 

The advantage of using the volume entropy is that it 
provides a good mechanical analogue of thermodynamic 
entropy even for small system, thus accounting properly 
for the finite-size effects. As the development of technol- 
ogy is allowing experimentalists to probe the thermody- 
namic behavior of smaller and smallers systems, this is 
becoming an increasingly important task. The main lim- 
itation of the present approach is that it is restricted to 
classical statistical mechanics, thus it does not account 
for quantum-mechanical phenomena. 
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